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Fully quantum treatment explicitly presents the high harmonic generation 

as a three-step process: (i) above threshold ionization (ATI) is followed by 

PLi| (ii) electron propagation in a laser-dressed continuum. Subsequently (iii) 

stimulated (or laser assisted) recombination brings the electron back into the 

initial state with emission of a high-energy photon. Contributions of all ATI 

channels add up coherently. All three stages of the process are described by 

^ \ simple, mostly analytical expressions that allow a detailed physical interpre- 

^^ I tation. A very good quantitative agreement with the previous calculations 

^' on the harmonic generation by H~ ion is demonstrated, thus supplementing 

,J^ ■ the conceptual significance of the theory with its practical efficiency. The 

, ^ . virtue of the present scheme is further supported by a good accord between 

the calculations in length and velocity gauges for the high-energy photon. 
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I. INTRODUCTION 

Under the influence of an intensive electromagnetic field an atom can emit electrons and 
photons. The number of photons absorbed from the field in the first process generally can 
exceed the minimum necessary for ionization resulting in distribution of the photoelectrons 
over the above threshold ionization (ATI) channels. The photon production manifests itself 
as the harmonics generation (HG) for the incident monochromatic laser radiation. Both ATI 
and HG are capable of populating the channels with remarkably high energy, as has recently 
been registered in experiments (see, e.g., Ref. [|l|-§]) and tackled by the theory ||6|-p3| (the 
list of references is unavoidably incomplete, for more bibliography see the reviews P^ , p5| ). 

An idea that the two processes referred above are interrelated was articulated long ago. 
Since in the HG process an active electron ends up in the initial bound state, it is appealing 
to represent it as ionization followed by recombination. This mechanism presumes a strong 
interaction between the emitted electron and the core that is omitted in the standard Keldysh 
p6| model of multiphoton ionization. The importance of this interaction was first pointed 
out by Kuchiev |2^ , who predicted several phenomena for which the electron-core interaction 



plays a crucial role. The related mechanism was named "atomic antenna". 

Specifically for HG, the simple relation between this process and ATI was suggested 
by Eberly et al ^ but proved to be non-realistic, see below. The hybrid classical-quantum 
model due to Corkum [§ (see also the paper by Kulander et al |jlO|]) casts HG as a three-step 



process: tunneling ionization and subsequent propagation in the continuum is completed by 
recombination. This intuitive model has influenced many research in experiment and theory. 
The simplicity of the model is due to some drastic presumptions. Usually it is emphasized 
that the intermediate electron propagation in the laser field is described by the Corkum |Q] 
model classically. Probably less attention is paid to the fact that neither the tunneling ion- 
ization through the time-dependent barrier, nor the laser-stimulated recombination receive 
a genuine quantum treatment as well. Being successfully applied to the comparison with 
some experimental data, the model resorts to such a loosely defined free parameter as the 
'transverse spread of the electron wave function'. From the conceptual side the Corkum 
model does not appeal to ATI process just because the discrete ATI channels do not appear 
within the classical framework. The subsequent theoretical developments were based on 
more sophisticated approaches and led to important advancements |TT|-|T^|K| , but appar- 
ently abandoned a perspective to establish a quantitative relation between ATI and HG. 
The ATI characteristics merely do not emerge in the papers devoted to HG theory, with 



few exceptions ||T^,p!6[]. For instance, Ref. |20[ establishes what the authors think to be 'the 



most general formal relation between ionization and high harmonic generation' that contains 
functional derivative of 'the ground state persistence amplitude' Z but not ATI amplitude. 
As the authors recognize, this relation is 'of limited practical use' since it requires knowledge 
of Z for arbitrary electric field F(t). 

Theoretically HG is governed by the Fourier components in the field-dresses wave func- 
tion. The proper description of the high-order components is a challenging task for the 
theory. Although some important features could be understood in classical calculations, 
the significance of fully quantum quantitative theory could not be overestimated. The 
computer-intensive numerical studies achieved substantial progress [0,||, but the available 
computer facilities often limit them to the one-dimensional models |6| j2T|j2^ . They could 



be favorably complemented by analytical studies capable of providing an important physi- 
cal insight. However, with the two different analytic quantum approaches being developed 
recently [0,|I^, the problem could not yet be considered as closed. For instance, as far 



as we know, these two approaches have never been applied to the same system in order to 
provide quantitative comparison of their results, albeit the differences in the formulation 



were indicated by Becker et al []15| . The only exception known to us is provided by Ref . pO| 
where in Fig. 6 some comparison is presented, albeit not for such a basic characteristics as 
HG rates, but for harmonic ellipticity and offset angle that reflect rather subtle peculiar- 
ities of HG process. For the offset angle a severe disagreement between the two theories 
is demonstrated. The lack of comparison looks quite typical to the current state of theory 
when the attention is mostly directed towards the qualitative aspects of the problem such 
as description of the pattern of emitted harmonic spectrum. 

The principal objective of the present study is to derive a fully quantum formulation for 
the HG amplitude in terms of the ATI amplitude and the amplitude of electron laser assisted 
recombination (LAR) in the laser field. Importantly, all the amplitudes are physical, i.e. 
no off-energy-shell entities appear. This circumstance adds to the conceptual appeal of the 
present theory its significance as a true working tool. In our approach, briefly described in 
the Letter to Editor |^, we base on the ideas outlined [^ and later developed in more 



detail [2£,^] by Kuchiev. We successfully test its efficiency by comparison with the bench- 
mark calculations by Becker et al |jl5| for HG by H^ ion. It should be emphasized that 



we test and achieve quantitative agreement for absolute values of HG rates, in contrast to 
many theoretical works which concentrate mostly on the qualitative issues. The quantita- 
tive character of our development is further illustrated by comparative calculations within 
the dipole-length and dipole-velocity gauges for the emitted high-energy photons: the good 
agreement testifies in favor of the method. In the broader perspective it should be empha- 
sized that our theoretical technique is directly applicable to other processes of current key 
interest, such as multiple ionization by laser radiation or enhanced population of high ATI 
channels due to the photoelectron rescattering on the atomic core. 

We start (Sec.||) with exposure of the general relations and outline our basic approxima- 
tions. In Section |T| we cast the harmonic generation amplitude as a sum over contributions 
coming from different ATI channels. The form and interpretation of this representation 
becomes particularly transparent when we resort to the adiabatic approximation (SecfTVp. 
The general theory is illustrated by the quantitative results in Section that is followed by 
the concluding discussion (Sec.|VI|). 

II. FORMULATION OF THE PROBLEM AND BASIC APPROXIMATIONS 

A. General relations 

The generation of a harmonic with the frequency Vt is governed by the Fourier transform 
of the dipole transition matrix element 

oo 



D{Vt)= dt exp{int) d{t) , (2.1) 

dit)= ldh^fir,tyde^iir,t) , 4 = e-r, (2.2) 



where \l/i and \I// are initial and final states of the atomic system dressed by the laser field 
(atomic units are used throughout the paper). 
We construct the initial field-dressed state \1/,- 



^,(r, t) = $,(r, t) + f dt' f dh' G(r, t; r', t') \/^(r', t') <l>,(r', t') 

developed out the initial field-free stationary state $„ 

<l>a(r,t) = ipa{r) exp{-iEat) , 



(2.3) 



(2.4) 
(2.5) 



using the retarded Green function G{r, t; r', t') which obeys equation 
d 



^g-^-Ha- VFit) 



G(r,t; r',t') = 6{t - t') 6{r - r') 
G{r,t;r',t')=0 (t < t') , 



it > t') 



(2.6) 



where Ha = 2P + ^a(r) is the Hamiltonian of an atomic system in the single active electron 
approximation, Va{r) is the interaction with the core, Vp(r,t) is the interaction with the 
laser field that generally includes the field-switching effects [below we presume that Vp(T) 
is real]. 

We consider harmonic generation when the atomic system ends up in the initial state. 
Other final states are also feasible, but till now this possibility have not been explored neither 
in the experiment, nor in the theory. In this case the final state tends to ^{t)a for t — > cx3 
being presented similarly to ( p.3|) as 



dt' J dh'<!>a{r',t'rVF{r',t')G{r',t'; r,t) , 
Employing (p.3| ) and ( |2.7| ) we transform d{t) to 

d{t) = {^a{t) I de I ^a{t)) + r dt' {^a{t) \ 4 G{t, t') W(t') | $a(t')) + 

+ / dt' ($„(t') I Vp{t') G{t', t) de I $„(t)) + 



(2.7) 



t 

+ / dt" / dt' {^a{t") I Vpit") G{t", t) de G{t, t') VF{t') I $,(t')) 

J t J —oo 



(2.8) 



where (. . . | ... | . . .) notation is employed to represent integration over the space variables. 



This formula could be compared with Eq. (2.15) in Ref. ||20| obtained under presumption 
\l/j = \l/j. Although the general structure is similar, the important difference lies in the 
range of temporal integration. In formula (|2.8|) the first term in the right hand side turns 
zero for the inversion-invariant (for instance, spherically symmetric) potentials. Physically 
the second term describes the process when the high harmonic photon is emitted after the 
interaction with the laser field, t > t', whereas the third term describes the 'time-reversed' 
event in which the radiation precedes the absorption of laser quanta. We denote the second 
and the third terms respectively as d'^{t) and d'it), 



rf+(t) = / dt'{<l>a{t) I deG{t,t')VF{t') I <l>,(t')) , (2.9a) 

J — CO 

/■oo 

d-{t)=J^ dt'{^a{t') I VF{t')G{t\t)de I $a(t)) . (2.9b) 



The last term in ( p.8| ) includes the effect usually referred to as continuum- continuum tran- 
sitions [llll , p!2| , p0| . In the present outlook it corresponds to the "mixed" picture when a part 
of laser photons is absorbed prior to emission of the harmonic followed by absorption of 
missing low-energy quanta. 

In this paper the continuum-continuum transitions are omitted as a rather standard 



approach, apparently assumed originally in Refs. |TT| , p!2| , albeit never scrutinized. There 
exists a simple physical reason why this term should not be important. The absorption of 
large number of low-frequency quanta happens when the active electron is well separated 



from an atom, see discussion of the role of large distances in p5|j36|l . In contrast, one should 
anticipate that emission of the high-energy quantum occurs when the electron is localized 
close to the atom. A transition of the electron from the outer region to the vicinity of 
the atom inevitably produces a suppression factor that describes the electron propagation. 
Later on it will be considered in detail, see the factor 1/ R in ( [4.14|) . For a natural sequence 



of events, when electron first absorbs A^ laser quanta and then emits the high-harmonic this 
suppression factor appears only once. For continuum-continuum transitions the electron has 
to go from the outer region into the vicinity of the atom thrice, that induces appearance of 
three suppression factors, thus substantially reducing the amplitude. Technically the explicit 
calculation of the continuum-continuum contribution means substantially higher level of 
difficulty, since it contains two Green functions; it is not pursued in the present study. An 
alternative, albeit indirect justification for omission of continuum-continuum contribution 
could be seen in the fact that our calculations of HG rates within this approximation are 
in a good quantitative agreement (Sec.0) with the results obtained by Becker et al who 
apparently do not rely on it. One more argument in favor of this approximation stems from 
the fact that it provides a very good agreement between the results in length and velocity 
gauges, as detailed in Sec. 0. 

Further we presume that the monochromatic laser field 

yp(r,t) = r-F coscut (2.10) 

is switched on adiabatically at some remote time (F is the amplitude of the electric field 
strength in the laser wave). In this case one can easily check that d~{t) = d'^{—t). After drop- 
ping in the right hand side of (2]^) the last term, that describes the continuum-continuum 



transitions, we arrive to the formula 

d{t) = d+{t) + d+{-t) . (2.11) 



One should remember that in our calculations according to (2^) and (|2.7| ) we employed 



different wave functions \l/i and \l//. Unfortunately, there is a strong trend in the literature 
H JTT|JT^JT5|J18|j3T| -p4| (to cite only part of references) to presume that \l/j = ^/ that leads to 
replacement of the relation (|2.11|) by 

d{t) = d+{t) + [d+{t)y , (2.12) 



and, consequently, to the real- valued dipole momentum d{t). 

Contrary to this almost universal delusion the expression ( p.l2 ) is incorrect for descrip- 



tion of harmonic emission by a single atom and the dipole momentum d{t) is not real- 



valued, both in the exact formulation (|2.8|) and within the approximation neglecting the 
continuum-continuum transitions ( |2.11| ). The recent paper on the unified theory of har- 
monic generation by Becker et al ||20| describes two different approaches. The calculations 
of the dipole-moment expectation presume \1// = \l/j and lead to the real-valued d{t) [see 
Eqs. (2.15) and (2.16) in |2y]], whereas the S'-matrix approach accounts for the distinction 
between \l/j- and \l/j and provides a similar formula (2.30) but with a different range of tem- 
poral integration [cf. our discussion below Eq. ( p.8| )]. The mentioned formula (2.30) gives 
complex-valued results. We agree with the authors comment that the S'-matrix approach 
is to be employed when the harmonic emission of a single atom is considered; namely this 



process is the subject of our present study. However, Ref. |2^ as well as Refs. [Il3| , p!^ carry 
out all calculations for the dipole moment expectation presuming that just this quantity is 
required as a source term for the integration of the Maxwell equations when the emission 
by the medium is considered. We believe that the theory of collective emission should be 
ultimately based on the proper description of a single atom process, but further discussion 
of this issue is beyond the scope of the present paper. 

For the high harmonic generation the 'time- reversed' process described by d^{t) = 
d'^{—t), i.e. emission of high harmonics followed by absorption of a large number of the 
laser quanta, is strongly suppressed as compared with the "natural" sequence of events rep- 
resented by d'^it), i.e. when at first a large number of laser quanta is gained and subsequently 
one high-frequency photon is emitted. Therefore the further approximation of Eq.( |2.11| ) 

d{t) = d+{t) (2.13) 

should work well. Below (Sec.0) we demonstrate by explicit calculation of the Fourier 
components that indeed the term d~{t) gives negligible contribution. 

For monochromatic laser filed ( p.lOj ) d{t) is a periodic function of time with the period 
T = 27f/u! where u: is the laser frequency. Introducing the Fourier transform 



1 r^ 
dN = 

we see that 



dN = - I dtexp{iNujt)d{t) (2.14) 

T Jo 



D{n) = 271 J2 S{^ - Nuj) dN . (2.15) 



N 



This stationary picture does not account for the depletion of the initial state by the laser- 
induced transitions. To describe this effect one could employ for \l/j and \l// the quasi-energy 
states with complex-valued quasienergies; however we do not resort to this complication 
below (another method to account for the depletion effects was suggested by Lewenstein et 

aim)- 

Thus the problem under consideration is reduced to calculation of the finite Fourier 
transform dN (2.14). One readily notices that 



dt, = d-_N, (2.16) 



and therefore 

whereas ( p.l2| ) leads to the distinct result 

d^ = dj^ + [d_j^ 



(2.17) 
(2.18) 



N 



< 



"at 



If 



As discussed above, for high harmonic generation we anticipate that dt_ 

the term d'^^ is negligible, then the formulae (|2.17|) and (|2.18| ) agree. The quantitative 

assessment for the quality of this approximation could be found in Sec.0. 



B. Keldysh-type approximation 

Our basic approximation is to neglect the effect of the atomic core potential Va in the 
time-propagator G{t,t'). A similar assumption underlies the Keldysh |2^ model, whose 
recent adiahatic modification p5|-p^] gives very reliable quantitative results for photode- 
tachment. A useful extension of the Keldysh model accounts for the Coulomb electron-core 



interaction p9| , |40| . The Green function within this approximation is straightforwardly rep- 
resented via the standard Volkov wave functions $p(r,t) [6{x) = 1, (a; > 0); 6{x) = 0, 
(x < 0)]: 



G(r, t; r', f) = -^9{t - t') J -^, $q(r, t) %{r', t') 

Explicit expression for the Volkov functions is conveniently cast as 
'^p(r,t) = Xp(r,^) exp (-iE^t) , 



Xp(r,t) = exp<^2 



(p + k,)r - f^ {E^{r) -E^)dr+^\ 



where the factor Xp(r, t) is time-periodic with the period T = 211/00 

F . 

kt = — sm ut , 

u 



E, 



2 



(2.19) 

(2.20) 
(2.21) 

(2.22) 
(2.23) 
(2.24) 



Due to the property ( p.l6|) it is sufficient to restrict subsequent analysis to the d'lj com- 
ponent. Using ( p.9a| ), ( |2.14| ) and ( p.l9| ) one can rewrite Eq.( |2.14] ) as 

dN = -fJdt I dt' J -^, i^ait) I e^^*4 I %it)){^^it')\VFit')\^ait')) = 

-00 ^ ' 



t 



- - y rft / dt' 





(i^q 



(</'a|4Uq(t))(Xq(t')|^F(0l¥'a)X 



X exp 



(277)3 
(Ea-E^ + n)t + l (E^ - Ea) t' 



(2.25) 



with Q = Nuo. In the latter representation the phase factors with the phases hnear in t and 
t' are exphcitly singled out; the remaining factor in the integrand is T-periodic both in t 
and t' . 



III. GREEN FUNCTION REPRESENTATION 

A. Factorization technique 

In the present section we use the theoretical technique elaborated by Kuchiev p9| , |30[l . It 
allows us to transform the right hand side of ( |2.25| ) to the form convenient for the analysis. 
This transformation amounts to some special representation of the time-dependent Green 
function, since ( |2.25| ) can be considered as its generalized matrix element. 

Generally speaking the correct description of the high Fourier components d^ represents 
a formidable theoretical task. Its numerical implementation via solving the non-stationary 
Schrodinger equation requires both a supercomputer and exceptional effort. In representa- 
tion ( p.25| ) the difficulty lies in the strong variation of the integrand as a function of the 
time variables t, t' . The application of the asymptotic technique is hindered by the fact 
that in expression (|2.25|) the integration variables are not independent: namely, the limit 
of t'-integration depends on t. The crucial simplification gained by using the factorization 
technique p9[ allows us to disentangle the integration variables at a price of introducing an 



extra summation. Very importantly, this summation is physically meaningful as it corre- 
sponds to the contributions of different ATI channels. The integration over the intermediate 
momenta q [coming from (|2.19|) ] is carried out in closed form. It should be emphasized that 
we do not use the so called pole approximation |^ applied recently by Faisal and Becker 



in their model of non-sequential double ionization by laser field pi| , ^ . 

In order to implement this program we transform the integral over t' using the identity 

/ dt'eM^Et') fit') = -^ Y. T^ I dt'fit') exp {z [{E - muj)t + mut']} , 

(3.1) 

that is valid for any periodic function f{t) = f{t + T). The identity can be easily derived 
with the help of the Fourier expansion f{t) = J2mfmQ^p{—irnujt). Employing it we rewrite 
Eq. ( p.25| ) in the form of the series 

djy = 2^dp^^ , (3.2) 

m 

where 

^"^ T^ Jo Jo J (27r)3 E^ - Ea - muj - iO ^ ^ 

X exp [i{Q — muj)t + imojt'] . 

The next step is to carry out the integration over q. To this end we note that the wave 
function Xq(r, t) [see Eq.(2.21)], depends on q in a very simple way. Namely, it is an exponent 



8 



of a linear form of q: Xq(r, t) = exp {iq [r + (F/u;^) cosut] + a}, where a is q-independent 
phase. The other source of q dependence in the integrand is its denominator that contains 
term |g^. The calculation of three-dimensional integrals of this type is standard. The 
integration over |q| might be carried out as a contour integration in the complex plane in 
order to specify the following validity conditions: 

ReR > , (3.4) 

ImKrn > , (3.5) 

where 

K^ = .2lmuj-^ + Eaj (3.6) 



is the photoelectron momentum after absorption of m laser quanta (see also Sec . [III B| ), 
R = vR? is the function of all the variables of integration 

R = R(r, r'; t, t') = I 'krdr — r' + r = ^r ( cos ujt — cos ut ] + r — r' . (3.7) 

Jt UJ^ ^ ' 

Note that below we shift the integration to the complex-valued time t', and hence complex- 
valued R. After integration over q the expression for (i^„ reads 



dtfm = -7^J dtj dt' J d^r Jdh'ipa{r){e ■r){r'-F COS oof) ipair') 



X 



2tiR 



exp < I 



KmR + i^- muj)t + mut' + k^r - k^/r' + [ dr i -k^ 

Jt' V 2 



4cj2 



(3.8) 



Eq.( |3.8| ) is convenient for evaluating the parameters governing the process. First of 
all note that integration over the variables r, r' is localized in the vicinity of the atom. 
The characteristic atomic dimensions should be compared with the amplitude F/uj"^ of the 
electron wiggling in the laser field (the latter motion being described by / hrdr). This 
amplitude becomes large even for quite moderate electric fields 

4 » 1 • (3-9) 

Note that this inequality may be satisfied both for large (7 > 1) as well as small (7 < 1) 



values of the Keldysh adiabaticity parameter 7 = ujJ2\Ea\/F. Therefore, for the fields 
satisfying ( p.9| ) we can assume that 

-^:>r,r'. 3.10 

Actually the applicability of the resulting approximation is even broader than outlined above 
since in fact relation ( |3.9| ) has to be replaced by its more accurate version 

(-eF) 

— I cosut — cosiut' 1^ 1 . (3-11) 



The main contribution corresponds to the complex values of time when | cosujt'\ ^ 1, as 
elaborated in the adiabatic approximation of Sec.|^ . 

These observations allow us to simplify R by neglecting r and r' in the pre-exponential 
factor of the integrand in ( p. 81) and retaining the first order terms in the phase: 



R^Rn + 



Rq • (r - rQ 
Ro 



Ro — Roi^, t ) — V R-o ! 



R-O — Ro(^;^ ) 



UJ^ 



(cos tut — COS tut') 



As a result the Green function in ( |3.8| ) becomes simpler in the 'wave zone 



— exp {iKrnR) ~ ^- exp {i [KmRo + K^ ■ (r - r')]} , 



K, 



F 



a = ±l . 



(3.12) 

(3.13) 
(3.14) 

(3.15) 
(3.16) 



Substituting Eq.( |3.15|) in (|3^ ) and using Eq. (|3.6|) to rewrite the exponent we find 
4m = - E ^2£dt£^ dt'Jdh Jdh'ifair) (er) (r'Fcosout') ^a(r') x 
— exp {i [KmRoit, t') + {VL — muj)t + mujt' + (K^ + k^) r— 



i?o(t,t 



Auo^ 



(3.17) 



Here summation over a = ±1 indicates that one has to take K^ parallel or antiparallel to 
F depending on the sign of (cos ut — cos tut') in order to satisfy the convergence condition 
( p.4|) (we postpone till the next section the more detailed discussion of this issue). Note that 
the denominator Rq in ( p.l7| ) can turn zero thus challenging convergence of integrals. This 
problem is circumvented if one presumes that the integration over t' is shifted from the real 
axis into the upper half-plane of complex-valued time. The integration contour C will be 
specified more exactly below. 

Now we return to the time-dependent bound state wave function $Q(r,t) and Volkov 
states $p(r, t) bearing in mind relation ( p.6|) . This gives us an appealing representation: 






E 



1 



27rT2 



dt / dt' 



Roit,t') 



X 



0(C) 

X ($„(t) I eMmde I $K^(t))(<l>K„(t') I VV(t') I $.(t')) . 



(3.18) 



B. Above Threshold Ionization and Laser Assisted Recombination 



An interpretation of formula ( |3.18| ) is based on the fact that the integrand is a product 
of physically meaningful factors. In order to recognize the first of them one should recall 



10 



that the amphtude of m-photon detachment of electron from the initial state $a within the 
Keldysh |2^ approximation is given by 

1 ^ 
Amip) = 7^jdt{%{t) I VF{t) I ^a{t)) . (3.19) 



Generally, the number m of photons absorbed is larger than a minimum necessary for the 
electron detachment, therefore the process corresponds to the Above Threshold Ionization 
(ATI). In the right hand side of (|3.19| ) the index m is implicit. It enters via the absolute value 



of the final electron momentum p which is subject to the energy conservation constraint 

hp' = mu-^ + Ea. (3.20) 

Note that —Ea = \n'^ > is the electron binding energy in the initial state; F"^ /[Aoj'^) is the 
electron quiver energy in the laser field. Comparing (|3.20|) with (|3.6| ) we confirm that Km 
is exactly the physical translational electron momentum in the m-th ATI channel. 

The second relevant process is that of laser assisted recombination (LAR) when the con- 
tinuum electron with the momentum p collides with an atom in the laser field. The electron 
absorbs some extra laser quanta and goes to the bound state $a emitting single photon of 
frequency Vt. In the Keldysh-type approximation this LAR process has the amplitude 

1 ^ 
CNn.{v) = -^^1^^^ ("^-W I exp(zfit) de \ %{t)) . (3.21) 



It could be called also laser-induced recombination. The allowed values of the high-energy 
photon frequency Q are 

1 F^ 

n = -p^ + -— + muj (3.22) 

2 4u;2 

with integer m. 

One readily notices that the integrand in ( |3.18|) bears a striking resemblance to the 



product of the integrands in ( |3.19| ) and ( ^.21| ). This allows us to interpret the amplitude 



( p.l8| ) as describing the two-step transition: at first the electron goes to the laser-field 
dressed continuum (Volkov) state ^k^ after absorption of m photons; at the second step it 
returns to the initial state emitting single photon with the frequency Q = Nuj. Note that 
the intermediate momentum p = K^ is restricted not only in magnitude (p.6|) but also in 



direction: according to Eq.( 3.16 ) the vector K^ is parallel or anti-parallel to the electric 
vector F in the laser wave. The physical implications of this circumstance are discussed in 
the next section. 

There is an extra factor 1/R{t, t') in the integrand of (|3.18|) which prevents complete 
separation of integrations in t and t'. In fact its presence is physically well understandable. 
Indeed, the definition of R{t, t') ([3171 ) shows that classically it is the distance between the 
electron positions at the moments t and t' with account for the electron wiggling in the laser 
field. 1/R could be named an expansion factor since in the absence of the laser field {F -^ 0) 
it describes conventional decrease of the amplitude in a spherical wave as it expands in 3D 
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space (~ l/|r — r'|). When the laser field is operative, the form of the expansion factor 
is drastically modified according to the approximation ( 3.14|) . Hence the interpretation of 
the expression ( |3.18| ) is that the electron first is transferred to the m-th ATI channel, then 
propagates in space under the infiuence of the laser wave and finally recombines to the initial 
state emitting the photon with the frequency Q. The contribution of each path is labeled 
by the number of virtually absorbed photons m. These contributions add up coherently as 
shown by Eg. ( ^.21 ). For the actual ATI process the momentum Km should be real. The 
summation in ( |3.2| ) includes also the virtual processes with the imaginary values of K^, but 
their contribution is anticipated to be small. 

It is worthwhile now to outline equivalent but more convenient representation of ATI and 
LAR amplitudes. Bearing in mind the character of r-dependence of Volkov states (|2.20|) - 
( p.21| ) it is easy to see that the space integration in the formulae ( |3.19| ) or ( |3.21| ) essentially 
reduces to the Fourier transformation. For instance, as shown in much detail by Gribakin 
and Kuchiev [^, the ATI amplitude ( |3.19| ) is exactly presented as 



A. 



[P) 



1 
'2T 



dt 



(P 



k,)' + ^' 



(p + kt)exp[i5(t)] 



(3.23) 



where S{t) is the classical action 

Sit) = ^fdT{p + Kf-Eat 



(q) is the Fourier transform of 0a(r): 



■a(q) = / dh exp(-zqr)0a(r) 



Similarly, for the LAR amplitude one obtains 



1 
^Nrn = -^^ jdt exp {i [fit - S{t)]} 0(^) (-K™ - k, 



(3.24) 



(3.25) 



(3.26) 



with 



(e), 



q) 



fe-V. 



iqj 



(3.27) 



Before concluding this section it is worthwhile to make an important observation. It is 
well known that the Keldysh approximation in the theory of multiphoton ionization is not 
gauge invariant. The detailed discussion of this issue can be found in Ref. |Q. Basing on 
physical grounds we use the dipole-length form ( p.lOj ) for the interaction of the electron with 
the laser field. As thoroughly discussed earlier [^^, this gauge stresses large separations 
of the active electron from the core where one-electron approximation is better justified and 
one can employ the asymptotic form for the initial-state wave function: 



AaT 



u~l 



exp(-Kr) Yim{r 



(r > 1/k), 



where k 



2\E, 



V 



ZJK^ Z is the charge of the atomic residual core 



V 



negative ion), / is the active electron orbital momentum in the initial state and r 



(3.28) 

for a 

= r/r is 



12 



the unit vector. The coefficients Aa are tabulated for many negative ions |^. The Fourier 
transform $a(q) ( P-25| ) is singular at q^ = 



with the asymptotic behavior for q —^ ±iK 



defined by the long-range asymptote ( |3.28| ) in the coordinate space 



{q)=47TAai±iyYirr,iq: 



(2/t)"r(z/ + i) 



2>+l 



(3.29) 



with (±1)' corresponding to q ^ ±in. 

At the present stage of our development the interpretation of ( p.l8| ) outlined above should 
be considered only as qualitative, since the two time integrations in Eg. ( p. IS] ) are not com- 
pletely separated due to the factor l/RQ{t,t'). The complete factorization is possible under 
additional approximation which actually is not restrictive possessing a broad applicability 
range. 

IV. ADIABATIC ANALYSIS 

A. Adiabatic approach to Above Threshold Ionization 



We start with reiterating the basics of the adiabatic approximation in multiphoton de- 
tachment theory ||35| , |36[| that allows one to carry out analytically integration in expression 
( p.l9| ) or ( p.23|) for the amplitude. It is presumed that the laser frequency u is small, i.e. 
that the number of absorbed photons m is large (the practical applicability of this approach 
proves to be very broad, since it gives reasonable results even for m = 2, see, for instance, 
Ref. [0). Then the integrand in ( |3.19[ ) or ( |3.23| ) contains large phase factor exp [^^(t)] and 
the integral may be evaluated using the saddle point method ||35|^8|. The positions of the 
saddle points in the complex t plane is defined by equation 



S'{t 



mfij 







(4.1) 



or, more explicitly. 



P + kt„^) + K^ = I p + — sin ut^^ J + K^ = . 



(4.2) 



Eqs. (|4.1| ) or ( [4.2|) are to be considered together with the energy conservation constraint 



Eq.( |3.20|) . Note that according to formula ( |3.29| ) the position of the saddle point coincides 
with the singularity of the bound-state wave function in the momentum space thus stressing 
the importance to describe correctly the long-range behavior of coordinate-space wave func- 
tion. The latter could be ensured much more easily than the proper description of the wave 
function inside the core region. Therefore the use of the adiabatic approximation and char- 
acterization of the bound wave function solely by its asymptotic behavior ( p.28 ) constitutes 
a self-consistent approach. 

The result of calculations of the amplitude ( |3.19| ) in the stationary phase approximation 



could be written as a modification of formula (25) in Ref. ||35|| : 
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^^"nP)=EAS(p), (4.3) 

^S(P) = -^ A. r(l + u/2) r'- .^ YUi>m,) I ^"P [^^(^"-^|] . (4.4) 

In the plane of the complex- valued time the saddle points tm/^ lie symmetrically with respect 
to the real axis. There are four saddle points in the interval < Retm^ < T, two of them 
lying in the upper half plane (Imtm^ > 0). The integration contour in the plane of complex 
time is shifted upwards. Therefore only two saddle points with Imt^^ > are operative 
being included into the summation in ( ^^ ) (/z = ±1); p^/^ is a unit vector in the direction 
of P + ki„^• 
B. Adiabatic approach to Harmonic Generation 



Although in the harmonic generation matrix element (|3.18|) the time integration variables 



t and t' are not fully separated, as indicated above, the large phase factor exp [iS{t')\ in t' 
integration is the same as in the multiphoton detachment case (|3.19| ). Therefore we can apply 
the saddle point approximation to carry out integration over t' . The integration contour is 
again shifted upwards in the complex t' plane; this is the contour C in the formulae (p.l7|) 



and ( |3.18|) . The saddle points are defined from the same Eq.(5]^). By solving it we find for 
the m-dependent saddle points t'^^\ 

sin ut'^^ = — {-Km + i^l^i) ■ 



cosut' = Wl - ( ^ ) {-Km + i^inY ilmcosut' > O) . (4.5) 



For each value of m there are two essential saddle points, i.e. these with Imt^^ > 0, labeled 
by subscript /^ = ±1. 

Thus in the stationary phase approximation for the t' integration we obtain from (p.l8|) 



dtimi^) = E E ^S(K«.) Bn^, , (4.6) 



T 



BNm, = -,--, fdt ($,(t) I exp(iOi) 4 I $K™(t)) . (4.7) 

Similarly to Eq.( |3.26D one can conveniently employ the Fourier transform which gives 

B,m, = -^jdt ^Xp{z[l]t-g(t)]} ~ _ 



rail 



In particular, for a negative ion {v = 0) with the active electron in an s state (/ = 0) we 
have from ( |3.29| ) (q = q/g) 

0a(q) = VA^Aa-^—^ , (4.9) 



2q_ 

' (g2 + fi:2 



^i^)(q) = -^{e. q) V^A^-—^^ , (4.10) 
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and ( |4.8| ) simplifies to 

T 



2Aa u^ f ,^ exp {t [nt - Sit)] } e ■ (K™ + k^) 



-Bat m/i = "^ ^~ ^T I dt 



v^T F J cos^t-cos^C^ [(K^ + ki)' + /€2 



(4.11) 



As noted above Eq.( |3.18| ), the summation over a = ±1 indicates that in order 
to satisfy the convergence condition ( |3.4| ) one has to choose a sign in RQ(t,t') = 
±(F/u;^) (coscut — coscut') that makes Ro(t,t') > 0. By considering the phase factor in 
( p.lTp , we see that the upper sign corresponds to K^ parallel to F whereas the lower sign 
leads to K^ antiparallel to F (without changing |Km|). This looks natural since there is no 
reason to prefer one of these directions. 



The right hand side of formula (|4.6| ) contains also summation over the saddle points 
t'^ (for the monochromatic laser field two such points are operative being labeled by the 



index /i = ±1). However, as argued by Kuchiev |^, actually only one saddle point (i.e 



one value of the label /x) contributes for each choice of a = ±1, i.e. for K^ parallel to 
F and for K^ antiparallel to F (as discussed below, this has a simple and clear physical 
interpretation). Consequently the double summation over a and /i is effectively replaced 
by a single summation. Moreover, in the latter sum both terms are equal. To see this one 
should remember that the preexponential factor l/RQ{t,t') changes sign depending on the 
value of cr = ±1. Hence the remaining single summation is equivalently replaced by the 
factor of 2. This allows us to finally rewrite ( p.2|) using (|4.6| ) as 



4 = 2 Y.Al:.ti^m)Br,^,,, (4.12) 

m 

where ^^^^^^(Km) and BNmfio are given respectively by ([4.4]) and ( [4.11| ) for negative ion 

(i/ = 0). Expression (|4.12| ) is to be calculated for the subscript fi corresponding to one of the 

saddle points t'^^=^^, for instance, that with the smaller value of Ret^^^^ (and ImtJ^^^ > 0). 

In formula ( [4.12|) the factor B^mfio describes jointly the 3D-wave expansion and LAR. 



These two effects could be further factorized using the approximation | coscut'^^J ^ | cos cut | 

Bntu^io = "^ CNm(K.m) , (4-13) 

4 = 25: Ai:%{Kj -^ CMm{K^) , (4.14) 

where 1/Rm^g = u'^/i^F cos t'^^^) is the laser-modified expansion factor in its simplest form 
and Cj\fm (|3.26|) is LAR amplitude. 

The adiabatic approximation could be further applied to carry out integration over t in 
BMmfj. ( [4.11|) or in CNm ( P-26D by the saddle point method. However below we do not pursue 
this objective and evaluate these integrals numerically. 

Eq. ( |4.12| ) and its simplified version ( [4.14| ) present the main result of this paper. These 



formulae implement the very simple picture of the HG process as consisting of three suc- 
cessive steps. First, the electron absorbs m laser photons. The amplitude of this event is 
Amu- In order to contribute to HG the photoelectron has to return to the parent atomic 

15 



core where LAR is solely possible. The amplitude of return is described by the expansion 
factor 1/R. It appears explicitly in ( 4.14 ), while in ( 4.12| ) it is incorporated in the definition 



of the amplitude B^mfio- The propagation of the electron describes the second phase of the 
event. At the third step the electron collides with the core absorbing N — m photons from 
the laser field and emitting the single high-frequency quantum Q = Nuj as it recombines to 
the bound state. This LAR process has the amplitude CNm ( |3.21|) . The summation over 



m in the total amplitude d^ ( [4. 121 ) takes into account interference of the transitions via 



different intermediate ATI channels. 

This appealing physical picture is supplemented by the very simple way to evaluate 
numerically all the quantities in ( [4.12| ). The amplitude of photoionization A^^^ is calculated 



via plain analytical formulae with the validity well testified before. The LAR process did 
not attract much attention in the literature and certainly deserves more study that we hope 
to present elsewhere. Here we emphasize only that since CMm^io as well as closely related 
amplitude Bj^rn^io are very similar in structure to the amplitude ^^^^q, one can hope that 
similar methods of evaluation also produce reliable results. 

The nontrivial point in the presented picture is the probability for the ATI electron 
to return to the core. Intuitively, one could anticipate that such a process is suppressed, 
because the most natural behavior for the electron would be simply to leave the atom. The 
proper description of the suppression plays substantial role in the theory. According to the 
physical image of the ATI process worked out in the adiabatic approach ||3^, after tunneling 
through the time-dependent barrier the ATI electron emerges from under the barrier at 
some point which is well separated from the core. As a result this point becomes the source 
of an expanding spherical wave. This occurs twice per each cycle of the laser field, at the 
two moments of time t' when the source-points lie up and down the field F from the core. 
The interference of the two spherical waves originating from the two different source-points 
results in non-trivial patterns in the angular ATI photoelectron distributions obtained from 
in agreement with the available theoretical and experimental data. The 



probability for the ATI electron to return to the core from the source-point is governed by the 
expansion factor 1/R and by the direction of propagation. At each of the moments t' only 
one of the two possible directions of K^, labeled in ( |3.18| ) by a = ±1, results in the electron 
eventually approaching the core. For the opposite direction of K^ the electron recedes from 
the core and does not come back to recombine. In other words, for each direction of K^ only 
one of the two saddle points t'^^ contributes to HG. Since both values of a give identical 
contributions, summation over a simply gives an extra factor of 2 in ([4.12|). 



C. Choice of the gauge 



Calculation of d^ according to the formulae (2^), ( p.l4| ), with exact wave function \E' 



might be equivalently carried out in various forms which correspond to different choice of 
gauge. As soon as the approximations are employed, the theory looses gauge invariance. 
Although the length gauge is known to be superior for the description of ATI within the 



adiabatic approximation [^, the situation is not that straightforward for the high-energy 
photon. In the derivation above we employed the dipole-length form for the operator d^ 
(p.2|). In the velocity- length form this operator is to be substituted according to the rule 



[p = — zVr is the electron momentum) 
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4= e-r^ -e-p (4.15) 

when calculation of (In is concerned (Q = Nu). It is easy to see that this substitution is 
equivalent to replacement of B^ ^ (|4.11| ) by 

Q^TfJ cos^t-cos^t;,^ [(K,„ + k,)V«:2] 
As compared with the expression (|4.11| ), the latter one differs by the extra factor 



" *^ — (4.17) 



2f] 
in the integrand. 



V. RESULTS OF CALCULATIONS 

Within framework of the present theory we calculate the rates of generating the iV-th 
harmonic radiation (c is the velocity of light) 

7^7V = ^^ I dN I' (5.1) 

introduced by Becker et al ||T5[ (and denoted by these authors as dR^ /dVL-Kj who in particular 
carried out calculations for the HG by H~ ion in the uj = 0.0043 laser field. Fig. |l| provides 
a complete comparison of these results with the results of our calculations. We employ the 
binding energy of H~ (k = 0.2354), but replace the true value Aa = 0.75 [^ by ^/2k = 0.686 



since this corresponds to the zero-range potential model used by Becker et al [|l5l. Our rates 
were multiplied by the extra factor N^, where N^, = 2 accounts for the presence of two active 
electrons in H~. For the real H~ ion the results shown in Fig.|| are to be scaled by a factor 

Fig. ID shows harmonic spectra for the laser field intensities / = 10^°, 2 ■ 10^*^, 5 • 
10^°, 10^^ W/cm^ that correspond respectively to the Keldysh adiabatic parameter values 
7 = LOKJF = 1.898, 1.342, 0.849, 0.600. Thus the most interesting region of transition 
from the multiphoton regime (7 ^ 1) to the tunneling mechanism (7 <^ 1) is covered. Our 
major results are shown by open circles in Fig. |1|. They are obtained using the expression 
( |4.12| ) for djy with the time-integration in B^mfio ( [1-11|) carried out numerically. Based on 



physical arguments, we extend the summation in ( [4.12 ) only over open ATI channels with 



the real values of K^- Generally the HG spectrum is known to consist of the initial rapid 
decrease, the plateau domain and the rapid fall-off region. The present theory is designed to 
describe the high HG but not the initial decrease, which in the case considered is noticeable 
only for one or two lowest harmonics. In the fall-of region, i.e. on the large- A^ side, our rates 
perfectly coincide with those obtained by Becker et al [jl5| (closed circles in Fig. |l|). The 
slight difference that could be hardly distinguished in the plot scale lies within uncertainty 
in retrieving data from the small-size plot published in Ref. |T^. Of course, our comparison 
is carried out in absolute scale without any fitting or normalization. 

17 



The deviations increase as N decreases, but within entire plateau region the agreement 
of rates averaged over structures remains good. Remarkably, the positions of numerous dips 
and peaks that exist in the plateau region are well reproduced by our calculations, although 
there exist some, generally not strong, discrepancies in their magnitudes. The structures 
has not yet received a universally accepted physical interpretation in the current literature 



with two tentative explanations being available. Becker et al 0] relate these structures 



to ATI channel closing whereas Lewenstein et al EUf suggest that it stems from quantum 



interferences between the contribution of different electron trajectories. Within the present 
framework we can say that the origin of these structures lies in the interplay of interfering 
contributions of various ATI paths, but their precise description would require an additional 
detailed analysis. 

Approximation ( |4.14| ) with numerical calculation of the integral in CNm ( p.26| ) (open 



squares in Fig. |T|) somewhat overestimates HG rate, but still retains the structure, though 
smoothed. It is worthwhile to emphasize two circumstances. First, to the best of our 
knowledge, we present here almost unique quantitative comparison of HG rates calculated 
within different theoretical methods, namely, by the present approach and by that developed 
by Becker et al. We are aware of only two cases when the quantitative comparison was carried 
out previously, one refers to some parameters of HG by elliptically polarized light in Ref. ||2^ 



[see also discussion in the Introduction), the other is comparison between the approach of 



Ref. 1 12 1 and results of direct numerical integration of the three-dimensional time-dependent 



Schrodinger approximation ||5T| . Second, the comparison is presented in the log-scale, as are 
the results reported by Becker et al, because it is appropriate both for the physics of the 
problem and for the current state of experiment. 



In the summation (^4.12 ) over ATI channels (i.e. over m) the coherence is very important, 



since large number of terms is comparable in modulus, but have rapidly varying phases. 
Many ATI channels contribute to HG for each N. This is in variance with the tentative 
conclusion by Eberly et al Q). The low ATI channels give appreciable contribution even 
for quite high harmonics. Only for the highest harmonics considered the contribution of low 
ATI channels becomes negligible. 

Although the length gauge is known to be superior for the description of ATI within the 



adiabatic approximation |^5[, the situation is not that straightforward for the high-energy 
photon. Therefore our calculations for the rates were carried out using both the length and 
velocity gauges. For large N the results obtained are very close, see Fig. ||; the dipole- velocity 
gauge producing the rates about 10% larger. For smaller N the difference increases and 
manifests rather irregular iV-dependence. However, even in the most unfavorable situation 
the ratio of the length-form to velocity-form results deviates from 1 not more than by 50%, 
which is rather reasonable, bearing in mind the multiphoton nature of the process. 

An evolution of the parameters of the individual A^-th harmonic radiation with variation 
of the laser field intensity / is presented in Fig. Rl for some particular values of A^. We show 

2 

both d^ (that is proportional to HG rate) and the phase ^n = argd^. The intensity 
dependence of the phase $Ar is known ||47|-|50|,^ to play an important role in description of 
the harmonic field propagation in the experimental conditions when the spectral and spatial 
coherence properties are substantial. Our calculations demonstrate rapid variation of the 
phase $Ar with the intensity J: as I increases by an order of magnitude the phase changes by 
about lOvr (for A^ = 9) or ISvr (for A^ = 11). The dependence of phase on I looks most simple 



in the fall-off domain where it is essentially linear, Fig. ^. As it is well known [p^p!5| , the 



HG rate manifests spikes (see Fig. ^, where the spikes are slightly smoothed due to finite 
step over / used in plotting) at the intensities /^ that correspond to the threshold of m-th 
ATI channels, that is closed for I > I^ due to ponderomotive potential F'^/{4uj'^)). 

As discussed in Sec.fl|, in the Fourier component dN ( |2.17| ) the contribution of rf^jy = '^jv 
is anticipated to be negligible as compared with that of d^. The reader have to recall that 
the physical reason for this is that d^ describes the natural sequence of events, when electron 
at first absorbs energy from the laser field and subsequently emits the high harmonic photon, 
whereas for d]^ the sequence is inverted. In order to illustrate how strong the preference is 

2 

we show in Fig. |^ the ratio djj/d^ . The ratio rapidly decreases with A^ being very small 
in all cases where the present theory applies. It becomes noticeable only for small A^ in case 
of intense laser field. 



VI. CONCLUSION 

The major result of this paper is a quantum mechanical description of the high-harmonic 
generation problem as a three-step process. In the first step the atomic electron absorbs 
several laser quanta and populates some ATI channel. Secondly, the electron propagates in 
the laser field back to the atomic core. The third step is the laser-induced recombination 
when the high harmonic photon is emitted. This mechanism is nothing more but interpre- 
tation of our principal formulae ( [4.12|) , ([4.14|) derived quantum mechanically with minimal 



approximations and without resort to any classical or intuitive arguments. The distinctive 
feature of our formulae ( [4.12| ), ( [4.14|) is that they include only genuine amplitudes of the 
three constituent processes, i.e. each amplitude describes true physical fully accomplished 
process, as HG, ATI or LAR. Here lies a crucial difference between our results and those of 
previous authors 11^ (see also Introduction) who discussed mathematical structure of the 
matrix elements or integrals considered in their approach in terms of three-step mechanism 
but failed to present HG amplitude via amplitudes of physically observable accomplished 
ATI and LAR processes. Even in a more broad context this is a rare and remarkable situa- 
tion when a complicated process is reduced essentially exactly to the sum over all possible 
paths with every path described as a sequence of real, fully accomplished physical processes. 
This conceptual simplicity arises due to multiphoton, adiabatic nature of HG process. 

With all three constituent amplitudes available from analytical formulae or simple nu- 
merical calculations our theory promises to be an efficient practical tool. This hope is 
supported by good agreement of our quantitative results for HG by H~ ion with the previ- 
ous calculations by Becker et al |jl5| in a wide range of laser intensities and frequencies of 



the emitted quantum. As the simplest example of possible future extensions we note only 
that our approach can be straightforwardly applied to negative ions with the outer electron 
having non-zero orbital momentum, such as halogen ions. These species could be easier 
accessible for the experimental studies. 

The structure of the formulae ( [4.12|) , ( |4.14|) is so simple that it is tempting to suggest 



that they will work well if one substitutes in them ATI or LAR amplitudes obtained in 
some approximations other than those used in the derivation above. For instance, if ATI 
amplitude more accurate than that given by the Keldysh approximation is available from 
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some theory, one can employ it for HG calculations via ( ^4.12|) , ([4.14|) . This gives a hope to 
relatively simply improve account for the electron rescattering effects in the HG theory. 

The distinctive feature of the present theory of HG is the use of the representation based 
on the discrete set of ATI channels. To the best of our knowledge this simple and apparently 
straightforward idea has not been exploited before. It would be underestimate to consider it 
merely as a detail of theoretical technique. Indeed, if one pursues the objective of the most 
direct and far-reaching quantum implementation of the three-step mechanism of HG, then 
the use of ATI channel representation becomes an unavoidable and crucial point. Otherwise, 
if the amplitudes of ATI do not appear in the theoretical scheme, HG cannot be properly 
described as a three-step process. 

In this paper we do not discuss qualitative features of HG spectra, such as extension of 
plateau domain etc. Detailed discussion of these issues could be found in other publications 
PJTTjJT^, in particular in the paper by Becker et al |]15| whose results for HG by H~ ion 



we reproduce closely within our theory. This implies that the analysis of the numerical 



results carried out by Becker et al |jT^ is fully applicable in our case. Concerning the 
mechanisms and physical interpretation of HG spectra, the present theory, hopefully, can 
add more to our understanding. However, the related analysis requires further theoretical 
developments which would overburden the present, already quite a long paper. We hope 
to present these subsequent developments elsewhere. In a broader perspective, one can 
expect that modifications of our approach could be applied to a variety of processes such as 
population of high ATI channels and multiple ionization of atoms. 

As a summary, the three-step mechanism of the harmonic generation is ultimately justi- 
fied. The theory is quantitatively reliable and easy to apply. It gives an important physical 
insight being a particular realization of the general atomic antenna mechanism. 
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FIGURES 

FIG. 1. Harmonic generation rates ( ^.ip (in sec~^) for H^ ion in the laser field with the 
frequency to = 0.0043 and various values of intensity / as indicated in the plots. Closed circles - 



results obtained by Becker et al |15|, open circles - present calculations in the dipole- length gauge 
using the expression ( 4.11| ) for B^m^i, open squares - same but with the simplified formula ( [4.14 ) 
for B^mfi- 



FIG. 2. Ratio of rates for HG by H ion in the laser field with the frequency to = 0.0043 
calculated within the present approach using the velocity and length gauges, (a) - for the rates 



calculated using the expression ( 4.11 ) for B^mfi'- circles - / = lO^''; squares - / = 2- 10^'^; triangles 
- / = 5 • lO^''; diamonds 
(PI ) for BNn.^.. 



I = 10 W/cm ; (b) - same as in (a) but with the simplified formula 



FIG. 3. Evolution of HG parameters with variation of the laser field intensity /. The harmonic 



intensity parameter M^v = 2 log 



10 



4 



and the reduced harmonic phase ^n/t^ with $7v = argd^ 
are shown for harmonic of various order A^. The bars with numbers m indicate the threshold 
intensities Im such that for I > 1^ the ?7i-th ATI channel is closed due to ponderomotive potential. 



FIG. 4. Parameter Qiv = 21og]^o ^n/^n ^^ ^ function of harmonic number N for HG by H 
ion in the laser field with the frequency w = 0.0043 and various intensities: circles - I = W^^; 
squares - / = 2 • lO^*'; triangles -1 = 5- 10^'^; diamonds - I = 10^^ W/cm^. 
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